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ABSTRACT 

This  report  describes  and  discusses  some  of  the  most  common  and  widely- 
used  methods  for  analyzing  travel  time  data  to  obtain  the  layer  and 
velocity  structures  in  marine  subbottom  sediments.  Two  methods  for 
analyzing  travel  time  data  of  arrivals  reflected  from  subbottom  layers 
to  determine  interval  velocities  and  layer  thicknesses  are  discussed 
and  shown  to  give  acceptable  resvats  when  applied  to  a range  of  differing 
■water~sediment  thickness  ratios.  A method  for  analyzing  inteirval 
velocities  to  yield  sound  speed  as  a function  of  depth  is  presented 
and  shown  to  be  capable  of  accvurately  determining  subbottom  velocity 
profiles  under  restrictive  conditions.  Two  methods  which  may  be  used 
in  a single  sediment  layer  configuration  are  also  presented;  both  of 
these  methods  are  based  on  the  effects  of  so\md  energy  which  has 
refracted  in  the  subbottom  and  are  capable  of  determining  the  initial 
sound  velocity  and  gradient  in  the  single  layer. 
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The  purpose  of  this  report  is  to  discuss  and  to  describe  different  ^ 

analytical  methods  for  determining  velocity  profiles  in  unconsolidated 
sediments  from  experimental  data.  This  report  is  by  no  meains  exhaustive; 
its  intent  has  been  to  assemble  in  one  place  those  methods  which  would 
be  most  useful  to  underwater  sound  propagation  modelers  and  best  suited 
to  the  types  of  data  that  they  will  usually  have  access  to.  ^ 

J 

There  are  different  ways  of  obtaining  information  about  the  velocity  \ 

structiire  of  the  subbottom.  In  shallow  water,  divers  with  probes  can  i 

take  samples  of  the  upper  few  meters  of  sediment  for  later  lab  emalysls.  \ 

A somewhat  deeper  probing  of  the  sediment  structure  is  possible  with  J 

sediment  coring  appeuratus  deployed  from  a research  vessel.  The  cores 
thus  obtained  may  be  analyzed  in  the  laboratory  or  the  velocity  profile 
of  the  sediment  may  be  measxired  in  situ  with  a device  such  as  the 
ARL  Compressional  Wave  Profilometer,^  which  can  be  attached  to  a corer. 

These  direct  methods  are  limited  by  the  depths  that  a corer  can  probe, 
or  roughly  30  meters.  The  only  other  direct  methods  involve  the  analysis, 
either  in  situ  or  laboratory,  of  drilling  cores.  This  method,  though  i 

very  useful  in  many  acoustic  studies,  relies  on  well  cores  which  are 
often  widely  spaced  and  must  be  supplemented  with  other  methods. 

The  methods  addressed  in  this  report  are  Indirect,  acoustical 
techniques,  which  more  often  than  not  provide  the  only  way  of  estimating 
the  velocity  structures  of  the  unconsolidated  sediments  to  the  depths 
required  by  underwater  propagation  modelers.  These  indirect  methods  ^ 

will  be  concerned  only  with  the  unconsolidated  sediments  though  some 
of  the  methods  can  be  and  are  used  in  seismic  profiling.  Selsmological 
methods  such  as  those  involving  lateral  and  shear  waves  axe  not  discussed 
in  this  report. 

i 
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This  report  is  primarily  intended  to  serve  the  needs  of  those 

working  in  the  field  of  underwater  sound  propagation  modeling  who  have 

little  knowledge  of  the  techniques  for  extracting  subbottom  velocity 

information  which  is  necessary  input  for  their  models.  An  example  of 

o 

such  a usage  is  the  ARL/UT  propagation  studies  using  BOTLOSS,  a 
program  that  requires  the  sediment  velocity  profile  as  input.  BOTLOSS 
calculates  a plane  wave  reflection  coefficient  for  a multilayered 
sediment  structure  overlying  a solid  subctrate  and  is  used  in  a ray 
theory  propagation  model  and  in  various  sensitivity  studies.  It  is 
applications  such  as  these  that  the  report  addresses. 

The  different  analysis  techniques  discussed  in  this  report  may  be 
divided  into  three  categories  which  are  treated  in  Chapters  II,  III, 
and  TV. 


A.  Reflected  Arrival  Methods 


Chapter  II  treats  the  reflected  arrival  methods,  which  are  of  two 

types,  the  T (X)  methods  and  the  ray  parameter  methods.  There  are 
2 

three  T (X)  methods:  the  reduced  time-reduced  distance  method,  the 

3 U ••  5 

small  range  method,  ’ and  Durbaum's  Method. 


All  of  the  T (X)  methods  are  derived  under  the  assumption  that  the 

water  and  sediment  layers  in  question  are  isovelocity.  This  assvimption 

2 

has  the  effect  of  restricting  the  T (X)  methods  to  small  water  layer- 
sediment  layer  thickness  ratios  (less  than  15:1).  Each  of  these  methods, 
though  they  approach  the  solutions  in  different  ways,  relies  on  the  fact 
that  the  relation  between  the  travel  time  of  a ray  reflecting  in  a single 
isovelocity  layer  €uid  the  range  attained  while  traveling  through  that 
layer  is  a hyperbolic  one  (T  =A+BX  ).  “ 


The  T (X)  methods  develop  techniques  for  determining  the  velocity 
of  sound  in  a layer  and  the  thickness  of  the  Is^rer  by  analyzing  travel 


time  curves  of  ray  paths  that  have  reflected  off  the  bottom  of  the  layer 
in  question.  (A  travel  time  curve  is  simply  a curve  depicting  the  rela- 
tion between  the  range  and  the  travel  time  required  to  attain  that  range 
of  rays  traveling  between  source  and  receiver  positions  at  fixed  depths.) 
When  the  layer  in  question  is  underlying  other  layers  it  is  necessary  to 
subtract  the  effect  of  the  overlying  layers  so  that  the  hyperbolic 
relation  between  T and  X in  a single  layer  may  be  used  to  yield  the 
layer  velocity. 

In  the  text  of  Chapter  II  the  reduced  time,  reduced  distance  method 
and  the  small  range  method  are  illustrated  with  examples.  Both  of  these 
methods  are  used  ^o  find  the  layer  velocity  in  a 200  m thick  sediment 
layer  underlying  a 500  m thick  water  layer  (thickness  ratio  2.5:1). 

The  reduced  time,  reduced  distance  method  determined  a layer  velocity 
of  1539  m/sec  and  the  small  range  method  yielded  1575  m/sec.  The  standard 
value  calculated  from  the  results  of  Appendix  A was  156I  m/sec  which 
yielded  percent  differences  of  1.4  and  O.89,  respectively. 

The  ray  parameter  method,^  which  is  the  second  of  the  reflected 
arrival  methods,  is  fundamentally  different  from  the  T (X)  methods, 
though  it  too  analyzes  travel  time  data  of  ray  paths  which  have  reflected 
off  subbottom  reflecting  horizons.  The  ray  parameter  method  is  unique 
in  that  water  or  sediment  layers  overlying  the  layer  to  be  analyzed  are 
irrelevant  to  the  solution  so  that  the  velocity  structure  above  the 
layer  does  not  need  to  be  dealt  with.  Therefore  the  ray  parameter 
methods  can  be  used  for  much  larger  water  layer-sediment  layer  thickness 
ratios  (up  to  500:1). 

The  ray  parameter  method  is  so  named  because  it  compares  the  ray 
paths  of  the  same  ray  parameter  which  have  reflected  off  the  top  and 
bottom  of  the  layer  in  question.  The  ray  parameter  is  defined  by 
p=cos0(z)/c(z)  and  is  calculated  frcxn  travel  time  data  by  evaluating 
the  derivative  of  a travel  time  cxurve,  i.e.,  p=dT/dx.  Two  different 
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types  of  ray  parameter  methods  are  discussed  in  this  report,  the  'normal 
method'  and  the  'thin  layer  approximation'.^  Both  methods  calculate  the 
soxind  velocity  and  thickness  of  a particiaar  layer  hy  matching  the  deriva- 
tives of  travel-time  curves.  The  thin  layer  approximation  incorporates 
a graphical  technique  which  alleviates  instabilities  associated  with  thin 
layers. 

An  example  of  the  thin  layer  approximation  is  considered  in  Chapter  II. 
This  method  is  applied  to  a single  sediment  layer  50  m thick  underlying 
a 5000  m water  layer  (thickness  ratio  100:1).  In  this  example  the  thin 
layer  approximation  yielded  a layer  velocity  of  I5U9  m/sec,  which  compares 
favorably  with  the  value  calculated  from  the  results  of  Appendix  A 
(1565  m/sec). 

B.  The  V(t)  Method 

Chapter  III  discusses  a method  whereby  velocity  measurements  obtained 
by  reflected  sirrival  methods  may  be  analyzed  to  construct  a velocity  profile 
throughout  the  sediment  layers.  As  stated  earlier,  the  reflected  arrival 
methods  determine  a layer  velocity  under  the  isovelocity  assumption,  a 
condition  which  is  rarely  met  in  the  physical  world.  For  this  reason, 
the  velocities  calculated  by  the  methods  of  Chapter  II  are  interpreted 
as  'interval  velocities'  which  are  assumed  to  be  the  time  averaged  sound 
velocity  in  a particular  layer  or  interval. 

T 8 

The  V(t)'’  method  is  a method  for  analyzing  interval  velocities 
to  obtain  velocity  gradients  and  layer  thicknesses  throughout  the  sediment 
layer  structure.  The  method  is  based  on  a velocity  function  which  is 
formed  by  fitting  the  layer  or  interval  velocities  to  a polynomial  in  t, 
the  one-way  vertical  travel  time  to  the  time  midpoint  of  the  various 
layers  on  the  sediment  structure.  Once  the  velocity  function  is  known, 
a gradient  function  and  a depth  function  which  specify  the  gradient  and 
depth  as  a function  of  t may  easily  be  determined. 
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The  effects  of  velocity  discontinuities  across  layer  boundaries  and 
uncertainties  in  the  interval  velocities  are  shown  to  be  significant  in 


some  situations  emd  tend  to  have  an  adverse  effect  on  the  confidence  one 


has  in  the  results.  In  particular,  velocity  discontinuities  yielded 


layer  gradients  in  error  by  an  average  of  30%  as  opposed  to  average 
differences  of  1%  when  velocity  discontinuities  were  absent. 


C.  The  Refracted  Arrival  Methods 

Q 

Chapter  IV  discusses  two  refracted  arrival  methods:  Dicus'  Method 

and  Hanna's  method. These  methods  differ  from  the  reflected  arrival 
methods  in  that  the  effects  of  shallow  refracted  arrivals  instead  of 
reflected  arrivals  are  analyzed;  also,  these  methods  assume  a gradient 
exists  in  the  sediment  layer.  It  is  expected  that  both  of  these  methods 
would  most  successfully  be  applied  to  thick,  single  sediment  layers. 

The  term  'refracted  arrival'  should  not  be  confused  with  arrivals  resulting 
from  lateral  or  interface  waves;  it  means  those  arrivals  that  are  returned 
from  the  sediment  layer  by  virtue  of  a turning  point  arising  because  of 
a velocity  gradient. 

Dicus ’ method  is  a graphical  method  which  deals  with  the  time 
differences  of  the  refracted  arrivals  and  the  arrivals  that  attain  the 
same  receiver  position  via  one  bottom  reflection.  The  method  is  derived 
assuming  a pseudolinear  sediment  velocity  profile  and  yields  the  gradient 
at  the  top  of  the  sediment  layer  and  the  ratio  of  velocities  at  the 
water-sediment  interface  as  output. 

As  an  example,  Dicus'  method  is  applied  to  a single  sediment  layer 
with  a constant  gradient  of  1 sec  having  a velocity  ratio  at  the  water 
sediment  interface  of  1.0333.  Dicus'  method  calculated  a velocity  ratio 
of  1.0358  and  a gradient  at  the  top  of  the  layer  of  0.99  sec”  , which 
correspond  to  differences  of  0.2k%  and  1%  respectively. 
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Finally,  a completely  different  method  due  to  Karma  is  considered. 
This  technique  yields  the  same  output  as  Dicus'  method  by  an  analysis 
of  the  nulls  in  propagation  loss  curves  which  arise  because  of  the  inter 
ference  produced  by  refracted  arrivals.  This  method  also  models  the 
sediment  with  a pseudol inear  velocity  profile. 


II.  REFLECTED  ARRIVAL  METHODS 
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This  chapter  discusses  the  most  common  and  widely  used  methods  of 
intervauL  velocity  analysis,  the  reflected  arrival  methods.  The  term 
reflected  arrival  methods  is  used  because  these  methods  deed,  with  arrivals 
that  have  reflected  off  the  bottom  or  subbottom  layers,  as  opposed  to 
methods  dealing  with  arrivals  that  have  refracted  through  the  bottom. 

Both  of  the  methods  to  be  discussed  amalyze  travel  time  curves  from  a 
layered  sovmd  medium.  A travel  time  curve  cam  be  constructed  for  each 
reflecting  surface,  and  this  travel  time  curve  is  simply  the  relationship 
between  the  horizontal  range  auid  the  time  required  for  a sound  ray  to 
attain  that  ramge  via  one  reflection  from  the  layer  in  question.  Though 
both  methods  approach  the  problem  in  different  ways,  each  extracts  the 
velocity  of  sound  in  the  layers  auid  the  thicknesses  of  the  layers  from 
travel  time  data. 

The  travel  time  data  to  which  these  methods  will  be  applied  are  all 
obtained  in  a similar  manner.  In  general,  a ship  deploying  TNT  charges 
or  towing  an  acoustic  source  travels  away  from  a configuration  of  hydro- 
phones which  are  either  physically  connected  to  another  ship  or  are 
attached  to  sonoboxiys.  As  the  source  ship  travels  away  from  the 
hydrophones,  the  charges  or  acoustic  sources  are  activated  at  given 
intervals.  Timing  circuitry  is  used  to  determine  the  travel  times  of 
the  arrivals  received  by  the  hydrophones,  and  this  travel  time  is  used 
in  conjunction  with  the  horizontal  distances  between  source  and  receiver 
to  construct  a travel  time  curve  (i.e.,  travel  time  versus  horizontal 
remge ) . 

When  charges  are  used  as  the  acoustic  source,  the  travel  time  data 
will  be  of  a discrete  nature.  When  sources  such  as  towed  sparkers  are 
used,  more  continuous  travel  time  curves  will  be  generated.  This 


L 
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continuous  data  is  usually  presented  in  the  form  of  amalog  paper  strip 
charts  that  €u?e  generated  on  hoard  ship.  These  strip  charts  give  a 
relatively  continuous  depiction  of  the  travel  times  of  the  various 
arrivals  versus  horizontal  distance,  or  ship  cruise  time. 


To  obtain  usable  mathematical  expressions  the  derivations 
discussed  in  this  chapter  axe  based  on  an  assvimption  of  isovelocity 
layers.  Unfortunately,  isovelocity  is  rarely  the  situation  in  the 
physical  world.  Most  often  this  method  is  used  to  find  the  sound 
speeds  in  the  subbottom  sediment  layers  of  a configuration  in  which 
the  top  layer  is  a relatively  thick  water  layer  and  the  underlying 
layers  are  thinner  sediment  layers  known  to  have  rather  significant 
sovind  speed  gradients.  The  sediment  layer  sound  speeds  will  be  inter- 
preted as  'layer  velocities*  or  'interval  velocities'.  These  are  the 
time  averaged  sound  speeds  in  the  layers  and  are  defined  i?y 
(see  Apx>endlx  A), 


C^(T)dT 
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where  AT  is  the  two-way  vertical  reflection  time  for  the  layer  and  C(T) 
is  the  soimd  speed  as  a function  of  time  for  the  layer. 

A.  The  T^(X)  Method 
2 

The  T (X)  method  is  one  in  which  the  square  of  the  arrival  times 
of  a reflected  arrived,  T,  are  fit  to  a polynomial  in  X,  the  horizontal 
reuige.  Usually  the  polynomial  in  X is  quadratic  but  more  advanced 
analysis  includes  terms  to  the  foiirth  power. 


2 11 
The  T (X)  method  is  based  on  analysis  done  by  Green  in 

If  one  considers  the  ray  path  through  an  Isovelocity  layer  as 


1938. 

depicted 
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in  Fig.  II. 1,  then  it  is  possible  to  write  the  following  Equation  for 
the  ray  path  length, 

+ (2H)^  . II. 1 

Since  the  layer  is  assumed  to  be  isovelocity,  the  above  equation  may  be 
divided  by  the  square  of  the  water  velocity,  C^,  to  yield 

T^(X)  = T^(0)  + -^X^  . II. 2 

"l 

2 2 
T (x)  is  the  reflection  travel  time  at  range  X and  T (O)  is  the 

vertical  reflection  time  at  the  origin.  From  the  form  of  Eq.  II. 2 it 

is  obvious  that  the  relation  between  T (X)  and  X is  a linear  one  with 

the  slope  being  C^”  and  the  intercept  being  T (O),  This  relationship 

is  exact  for  an  isovelocity  single  layer  configuration.  Green  first 

used  Eq.  II. 2 to  calculate  soiand  speeds  by  plotting  T (X)  and  X values 

obtained  from  seismic  measurements  and  measuring  the  slope  of  the 

resulting  graph.  When  the  configuration  is  changed  to  a more  realistic 

one  in  which  added  layers  with  differing  velocities  are  considered,  the 

2 2 

linear  relation  between  T (X)  and  X is  destroyed. 

To  illustrate  the  T (X)  method  and  its  different  versions,  consider 
the  config\aration  of  two  isovelocity  layers  depicted  in  Fig.  II. 2.  The 
ray  path  labeled  can  be  analyzed  using  Green's  method  since  its  path 
is  wholly  within  one  layer;  the  path  labeled  R^2  cannot  since  it  traverses 
two  layers.  To  handle  this  latter  situation,  Dix  (1955)  modified  Green's 
method. 

1.  The  Reduced  Time.  Reduced  Distance  Method 

To  calculate  the  velocity  in  layer  2,  information  from  paths 
like  R^2  must  be  analyzed.  As  mentioned  earlier,  for  paths  like  these 
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2 2 

there  is  no  longer  a linear  relationship  between  T (X)  and  X . However, 

in  many  cases  of  physical  interest,  this  linear  relationship  is  very  nearly 

3 U 

preserved.  The  method  developed  by  Dix  * uses  information  obtained  from 
paths  like  to  effectively  subtract  the  first  layer  from  the  problem  so 
that  Green's  method  can  thus  be  applied  to  the  'reduced  data. ' 

This  technique  is  as  follows.  As  stated  before,  the  T^(X)  vs  X^ 
relation  is  'quasilinear'  in  many  cases.  To  apply  the  reduced  time, 
reduced  distance  method  one  fits  the  T^(X)-X^  data  for  paths  like 
to  a straight  line.  When  this  is  done,  one  has  the  following  equation. 


2 1 ? 
T (X)  * B + X 


where  is  interpreted  as  a 'rms  average  velocity'  through  both  layers. 
In  doing  this  one  essentially  assumes  that  the  path  is  the  least  distance 
path  (dotted  line  in  Fig.  II-2)  rather  them  the  true  least  time  path 
specified  by  Fermat's  Principle. 

To  be  able  to  subtract  the  effect  of  the  first  layer,  the 
angle  6^^  is  needed.  In  most  cases  of  physical  interest,  this  angle  can 
be  approximated  very  well  by  0,  which  can  be  easily  calculated  frexn 
Eq.  II. 3 as  follows.  The  ray  pamameter  of  the  path  specified  by  Eq.  II. 3 
is  (see  Appendix  B) 


_ dT(X)  ^1 

■ dX  ' C, 


u 


Using  Eq.  II. 3»  one  obtains 


2T(X)  a 

dx  TT 


- <iT(X)  _ 
dX 


cos  e. 


T(X)V, 


or 


C X 

cos  0 « 


T(X)V. 


II. U 


Once  0^  is  known,  those  portions  of  the  time  and  horizontal  distance  of 


path  R^2  which  lie  within  layer  1 may  be  calculated  and  subtracted  from 
the  original  T and  X values  for  paths  is  done  in  the  following 

manner.  From  Fig.  II. 2,  the  contribution  to  X from  layer  1 is  AX,  where 


2Z, 


AX  = 


tan  0. 


II. 5 


The  contribution  to  the  time  resiilting  from  layer  1 is  At,  where 


2Z.. 


At  = 


sin  0^Cj^ 


II. 6 


£Uld 


2Z, 


sin  0. 


is  the  length  of  the  path  R^^  which  lira  in  layer  1.  When  AX  and  At 
are  known,  the  first  layer  may  be  subtracted  to  leave  the  reduced  time 
and  reduced  distance  values  and  X^^,  which  are  given  by 
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Tj^  = T - At 


= X - AX 


If  the  T and  X values  for  Paths  are  reduced,  the  problem  has 

essentially  been  transformed  to  a one  layer  problem.  This  being  the 

2 2 

case,  the  relation  between  and  should  be  a linear  one  given  by 


T ^ = (T  ) ^ 

R ' o^R  _ 2 ^ 

^2 


2 

where  (Tq)^  is  the  vertical  reflection  time  through  layer  2 and  is 

the  interval  velocity  in  layer  2.  Once  (T  )_  and  C„  are  known,  the 

On  d 

thickness  of  layer  2 may  be  determined  by  the  following  equation. 


This  method  is  best  Illustrated  by  considering  an  example. 

Figure  II. 3 depicts  the  T versus  X curves  for  a two-layer  situation  in 

which  the  first  layer  has  a thickness  and  velocity  of  500  m euid  1500  m/sec 

and  the  second  has  a thickness  of  200  m and  an  initial  sound  speed  of 

1365  m/sec  which  Increases  linearly  with  a gradient  of  2 sec”^.  These 

12 

curves  were  generated  from  a ray  tracing  program,  PROBE,  which  was 
developed  at  ARL/UT  by  Terry  Foreman. 


The  top  curve  of  Fig.  II. 3 results  from  R.„-type  paths,  whereas 

2 ■^‘^2 

the  bottom  one  resiats  from  R^^-type.  The  T versus  X plot  for  the  top 

p 00 

cxirve  is  shown  in  Fig.  II. U.  A least  squares  fit  yields  « 2.3l4  km  /sec  . 
If  Eqs.  II.  1*,  II. 5»  II. 6,  and  II. 7 are  used  along  with  X and  T values 
obtained  from  Fig.  II. 3.  then  Table  II. 1 may  be  computed.  The  Xp  and 
values  that  were  ccxnputed  are  then  fit  by  least  squares  (see  Fig.  II. 5) 
to  yield 


« 2.314  kmV»*c2 
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20  40  60 

X*  X 10~*  - km* 


J 

80 


FIGURE  n.4 
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TABLE  II. 1 


II. 1 


Tp  = 0.0661+2  + 0.1+2216 

From  Eqs.  II. 8,  9,  and  10  one  obtains 


„ (0,,066l.2£0. 1.2216)^^^  , „_„8 


1.539  km/ sec 


The  actual  values  of  AZ  and  are  200  m and  I56I  m/ sec , where  the 
interval  velocity  of  layer  2 was  calculated  from  Eq.  A. 6. 


This  reduced  time,  reduced  distance  method  can  be  extended  to 
the  case  where  there  sure  more  than  two  layers  by  a repeated  application 
of  the  method.  For  example,  if  three  layers  are  being  considered,  the 
velocities  in  the  first  two  layers  could  be  found  using  the  techniques 
developed  in  this  section.  The  knowledge  of  these  two  velocities  suid  the 
ray  parameter  would  sillow  the  first  two  layers  to  be  subtracted  off,  again 
leaving  a one-layer  problem  that  is  easily  haindled.  In  this  manner,  the 
method  may  be  extended  to  the  n-layer  case. 

2.  The  Small  Rsmge  Method 

Although  the  method  Just  discussed  gives  good  results  in  most 

situations,  it  is  relatively  hard  to  use  since  a large  n+mber  of  X suid  T 

3 U 

values  must  usually  be  reduced.  Dix  ’ developed  an  approximation  which 
is  valid  for  small  range  values  and  gives  good  results.  This  approxima- 
tion is  derived  as  follows. 

From  Fig.  II. 2,  the  following  geometrical  relations  may  be 
obtained  for  path 


I ? 

* 

♦ 

I 

I 

! 

i 


= 2Z^  cot  0^  = cot  6^ 


= 2(Z^~Z^)  cot  0^  = C^Tg  cot  0^ 


X = X^  + Xg 


where  X^  is  the  contribution  to  the  range  resulting  from  layer  i and 
is  the  vertical  two-way  travel  time  thickness  of  layer  i.  If  the 
range  values  are  small  compared  to  the  thicknesses  of  the  layers,  i.e.. 


X « + AZ  , 


then  the  cotangents  of  the  angles  may  be  approximated  with  the  cosines. 
This  substitution  yields 


Xi  + Xg  - cos  0^  + CgTg  cos  8^ 


- cos  0..  P cos  0^ 
^1  ^1  " ^2  "2 


But  by  Snell's  Law, 


cos  0^  cos  0g 


hence 


/ 2 2 \ ®1 

Xi.X2  = (Ci 


Substitution  for  cos  6^  from  Eq.  II. U yields 


or  as  X-*-0, 


T(X)V, 


£ 2 2 

where  V (O)  is  the  slope  of  the  T versus  k plot  for  R, „ paths  in  the 
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limit  of  small  range  values  and  T(0)*T^+T2- 

This  method  will  now  be  applied  to  the  same  example  considered 

in  the  preceding  section.  From  Fig.  II. U,  the  slope  of  the  versus  X^ 

c\irve  for  small  range  values  is  fo\and  to  be  V.^*2.3lU  km^/sec^.  The 

A 

velocity  In  the  first  layer  is  assumed  to  be  known  and  is  C^=1.5  km/sec. 
From  Fig.  II. 3 the  vertical  two-way  time  thicknesses  of  layers  1 and  2 
may  be  measured  to  give 


T^  = 0.6668  sec  , 


= 0.2577  sec 


Equation  11.12  may  now  be  solved  for  the  layer  velocity,  C, 


2 V^^(0)(T^+T2)  - 


= 2.3593  km^/sec^ 


Cg  = 1.575  km/sec 
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This  velocity  compares  favorably  with  the  results  of  the  previous  method 
(Cg=1.539  km/sec)  and  the  value  calculated  from  Eq.  A. 6 (C2*1.56l  km/sec) 
The  thickness  of  the  layer  may  be  computed  as  follows. 


This  resiilt  Is  also  In  good  agreement  with  the  true  value  of  0.200  km 


The  small  range  approximation  considerably  reduces  the  effort 
required  to  calctilate  interval  velocities  while  having  very  little  effect 
on  the  results.  This  method  is  very  well  suited  for  work  in  the  field 
since  the  results  can  be  obtained  quickly  using  graphical  techniques. 


This  method  may  also  be  used  in  the  n-leiyer  case.  By 
mathematical  induction,  Eq.  11.12  may  be  extended  to  the  n-layer  case 
to  give 


The  interval  velocity  of  the  nth  layer  is  therefore 


where  is  the  two-way  vertical  time  throiigh  the  ith  layer,  and  V ^ 
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is  the  reciprocal  slope  of  the  T versus  X curve  for  the  nth  reflector. 
3.  Durbaum's  Analysis 

2 

This  final  version  of  the  T (X)  method  is  a method  developed 
by  Durbaum^  in  195*»*  The  material  used  in  this  report,  however,  was 
tadcen  from  a paper  by  Clay  and  Rona^^  ( 1965 ) , because  it  provided  a 
good  summary  of  the  method  in  English. 

Durbaum's  method  is  a higher  order  method.  Instead  of  fitting 
travel  time  data  to  a quadratic  function  of  X,  the  data  is  fit  to  a 
fourth  order  polynomial  in  X,  i.e.. 


Tj^(X)  » Aq(J)  + A^(J)X  + + A^(J)X^  + A^(j)X*^ 


II. lU 


The  subscript  J refers  to  reflection  times  from  the  Jth  layer  and  T.(X) 

V 

refers  to  a two-way  reflection  time  rather  than  a two-way  time  thickness. 
The  6ulded  first,  third,  and  fourth  power  terms  in  Eq.  II. ll*  allow  layer 
configurations  of  a more  complex  nature  to  be  treated.  One  important 
case  that  can  be  treated  by  Eq.  II. lU  is  the  one  in  which  the  subbottom 
layers  are  inclined  at  a small  angle  with  respect  to  the  surface  leiyer 
(see  Fig.  II. 6).  The  coefficients  may  be  expressed  in  terms  of  the 
physical  euid  geometrical  parameters  of  the  layered  configuration  as 
follows , 
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CONFIGURATION  FOR  SLOPING  SUBBOTTOM  LAYERS 

(AFTER  CLAY  AND  RONA) 
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- ^ - ^)l 

J r , 

J , 

B,,  = 5^  C JT  (0)-T  ,(0)1 
y ^ m [ m m-1  J 


To  apply  this  version,  as  with  the  others,  it  is  necessary  to 
know  the  layer  velocities  and  vertical  travel  times  for  each  layer  above 
the  one  in  question.  When  these  are  known,  along  with  the  fitting 
coefficients  Aj^(j),  it  is  possible  to  calculate  the  layer  velocity  of 
the  layer  in  question  plus  the  angle  of  inclination  of  the  subbottom 
layers,  y*  (Notice  that  when  y=0,  the  A^(J)  and  A^Cj)  terms  are  zero.) 


If  the  fourth  power  term  is  neglected  and  If  an  angle  of  inclination  of 
zero  is  considered,  Durhaum's  method  reduces  to  the  small  range  method. 
In  practice,  Durhaum's  method  is  most  useful  when  sloping  bottom  layers 
are  present.  When  the  subbottom  is  not  sloping,  very  little  is  gained 
by  including  the  higher  power  terms.  For  more  details  concerning  the 
application  of  this  method,  consult  Clay  and  Rona^^  (1965). 
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U . Limitations  of  the  T (X)  Methods 


The  main  factor  which  causes  the  T (X)  methods  to  break  down 
is  the  violation  of  the  conditions  assumed  when  deriving  the  equations 
to  be  used.  The  most  apparent  condition  which  is  violated  is  the  iso- 
velocity assumption.  When  dealing  with  physical  situations,  the  isovelocity 
condition  is  almost  never  satisfied.  Hamilton^**  has  determined  that 
sediment  velocity  gradients  generally  range  between  O.5  euid  2.0  sec”^; 
hence  errors  are  introduced  because  the  paths  through  the  layers  are 
curved  rather  than  straight.  In  very  thick  nonhomogeneous  layers,  this 
can  render  the  T (X)  methods  useless  (Bryan,  197*t). 


As  a further  example,  consider  a configuration  like  the  one 
in  Fig.  II. 2,  where  the  top  layer  is  a relatively  thick  water  layer  and 
the  bottom  a much  thinner  sediment  layer.  Let  be  the  time  associated 
with  path  and  Tg  be  the  time  associated  with  isovelocity 

assumption  introduced  an  absolute  error  of  At  in  the  time  of  path 
then  the  relative  error  in  the  water  is  of  order  At/t^.  If  the  data 
concerning  the  second  layer  are  calcvilated,  the  relative  error  in  these 
data  is  of  order 


^2  ■ "^l 


"^l  "^2  " "^l 


Since  roughly  proportional  to  h^^/hg,  where  h^^  is  the  thickness 

of  layer  i,  the  absolute  error  of  the  data  in  the  bottom  is  magnified  by 


1 


F-! 


I 

i 
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the  ratio  of  the  thicknesses  of  the  two  layers.  In  most  cases  the  water 
layer  has  a fairly  slowly  changing  velocity  profile,  but  If  the  water 
layer  were  very  thick  and  the  sediment  layer  very  thin.  It  Is  conceivable 
that  the  absolute  error  of  the  reduced  data  could  be  of  the  same  magnitude 
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as  the  data  Itself.  These  are  the  conditions  under  which  the  T (X) 
methods  break  down. 

The  examples  considered  In  this  section  yielded  results  for 
layer  velocities  and  thicknesses  that  were  within  of  the  standard 

values.  These  results  are  quite  good,  but  It  should  be  remembered  that 
the  data  that  were  analyzed  were  obtained  from  computer  ray  tracing 
programs.  When  dealing  with  real  world  data,  the  quality  of  the  data 
will  usually  preclude  attaining  results  as  good  as  these. 

B.  Ray  Parameter  Method 

The  ray  parameter  method  (RPM)  developed  by  Bryan^  (l97l»),  uses  an 
alternative  approach  In  the  analysis  of  travel  time  data.  Rather  than 
compare  paths  of  the  same  shot  position,  the  RPM  compares  paths  of  the 
same  ray  parameter  p.  This  approach  yields  a solution  for  the  velocity 
In  subbottom  layers  which  Is  Independent  of  the  velocity  structure  above 
this  layer. 

1.  Derivation  of  the  Method^ 


A ray  path  Is  characterized  by  a constant  ray  parameter,  p, 
which  Is  defined  by  Snell's  Law  to  be 

cos  6(z) 

P = ~'cTzl  • 

The  angle  0 Is  the  angle  that  the  ray  makes  with  the  horizontal,  l.e., 
the  grazing  angle. 


i 

1 

j 

1 

I 

1 

I 

i 
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For  an  incremental  travel  time  dt,  the  following  equations 
may  be  written  (see  Fig.  II.7), 


dz  * c(z)dt  sin  9(z)  , 


dx  = c(z)dt  cos  9(z) 


If  Eq.  11.15  and  the  trignometric  relation 


cos^  9 + sin^  9=1 


are  used,  Eqs.  II. l6  may  be  written  as 


/ 22 

dz  = (1-p  c ) cdt  , 

2 / 2 

dx  = c pdt  = pc  1 1-p  cl  dz  , 


/ 2 2i 

where  cdt  = dz/yl-p  c j has  been  used  to  obtain  the  second  of  Eq.  11.17. 
The  travel  time  and  range  as  a function  of  p may  be  obtained  by  rearreuoging 
and  integrating  Eqs.  11.17  to  yield 


t(p)  = 2 


/ ./,  ^ 

o c yl  - P c 


.3/ 


11.18 


L 22 
'1  - p c 


To  transform  Eqs.  11.10  into  usable  forms,  consider  the 
configurations  pictured  in  Fig.  II. 8.  The  object  of  the  method  is  to 
measure  the  interval  velocity  of  the  layer  bounded  by  zs*z^  and  z»z^+Az. 
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dz  s c(z)  tin  9(z)dt 
dx  - c(z)  cot  9(z)dt 


FIGURE  n.7 

GEOMETRY  SHOWING  THE  RELATION  OF  dx  AND  dr  TO  dt 


ARL  ■ UT 
AS-74.U7* 
SRR  . DR 
II  • 2 - 7« 


29 


Il 


• -1 


The  two  paths  pictured  in  the  bottom  half  of  Fig.  II. 8 give  rise  to  the 
travel  time  curves  in  the  top  half.  Both  paths  have  the  same  ray  param- 
eter p.  This  is  a very  important  point  as  will  be  seen  shortly. 

To  calculate  the  layer  velocity  It  is  first  necessary  to 
evaluate  At  and  Ax  as  functions  of  p.  This  is  done  as  follows. 


The  fact  that  paths  with  the  same  ray  parameter  are  being  considered 
means  that  p in  Eqs.  II. 19  is  constant.  To  evaluate  Eqs.  11.19  only  one 
simplifying  assumption  need  be  made,  that  being  that  the  layer  bounded 
by  and  z^+Az  is  isovelocity.  If  this  assumption  is  made,  the  integrals 
in  Eq.  11.19  are  readily  evaluated  to  yield 


—1/  2 

At(p)  = 2C^  -^fl-p  Az 

/ p 

Ax(p)  = 2pC^^l-p  j Az  ; 11.20 

the  quantity  is  used  to  denote  the  interval  velocity  of  the  layer 
in  question. 
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Equation  11.20  may  now  be  solved  for  the  interval  velocity 
of  the  layer  to  give 


The  important  thing  about  the  RPM  is  that  the  velocity  structure 
above  z=z^  has  no  bearing  whatsoever,  a consequence  of  analyzing  paths 
having  the  same  ray  parameter.  Since  both  paths  have  the  same  ray  param- 
eter, their  paths  in  the  region  above  z*z^  are  completely  parallel  and 
hence  cancel  totally  when  the  integration  is  performed.  The  only  region 

where  isovelocity  was  assiuned  was  in  the  layer  to  be  analyzed.  Therefore 

2 

the  RPM  is  inherently  superior  to  the  T (X)  method  since  it  assumes 
isovelocity  in  all  layers. 


2.  The  Normal  Method  of  Solution 


To  apply  the  RPM  it  is  necessary  to  match  the  paths  of  two 
adjacent  reflectors  having  the  same  values  of  p.  The  data  to  be  analyzed 
will  be  travel  time  data  and  the  ray  parameter  will  vary  along  these 
curves.  The  parameter  of  a ray  arriving  at  a particular  range  is  the 
slope  of  the  travel  time  curve  at  that  reuige  (see  Appendix  B).  Hence, 
applying  the  RPM  will  involve  calculating,  numerically,  the  derivatives 
of  travel  time  curves. 


To  see  how  the  RPM  is  normally  applied,  consider  Fig.  II. 8. 
The  first  step  in  the  analysis  is  to  find  two  places  on  the  T-X  curves 
which  have  the  same  slope.  This  will  involve  calculating  nimierical 
first  derivatives  on  each  curve  and  matching  them.  Let  points  A and  B 
in  Fig.  II. 8 denote  two  such  points.  Next,  the  differences  in  times 
and  distance  between  the  points  can  be  determined.  'If  p,  Ax,  suid  At 
are  known,  the  interval  velocity  may  then  be  calc\ilated  by  the  first 
of  Eqs.  11.21.  When  C-  is  known,  Azmay  be  calculated  using  the  second 


of  Eqs.  11.21.  This  procedure  may  be  repeated  for  many  different  points 
and  the  results  averaged  to  give  the  layer  thickness  and  interval  velocity. 

Although  the  RPM  is  inherently  superior  to  the  T^(X)  method, 

it  is  not  without  its  shortcomings.  When  the  layer  to  be  analyzed  becomes 

thin  compared  to  the  overlying  layers,  the  RPM  becomes  unstable  as  does 
2 

the  T (X)  method.  This  instability  arises  for  different  reasons  though. 

When  a thin  layer  vuiderlying  a thick  layer  is  to  be  analyzed, 
the  T-X  curves  for  the  bounding  interfaces  will  be  practically  parallel 
over  the  entire  range  of  x.  This  means  that  points  A and  B in  Fig.  II-9 
will  have  almost  the  same  range  values,  and  this  introduces  an  instability 
because  Ax  is  so  small  that  the  error  associated  with  its  measurement 
is  of  the  same  order  as  Ax  itself.  Bryein  has  determined  that  these 
instabilities  arise  when  the  overlying  layers  are  roughly  50  times  the 
thickness  of  the  layer  to  be  measured.  Fortunately,  the  RPM  may  be 
modified  in  such  a manner  as  to  remove  the  instability  associated  with 
small  values  of  Ax(p). 

3.  The  Thin  Layer  Approximation 


The  thin  layer  approximation  (TLA)  may  be  derived  from 
Eqs.  11.20.  If  one  takes  the  ratio  of  At(p)  to  Ax(p),  then  the  following 
equation  is  obtained. 


= MeI 

_ 2 Ax(p) 


The  next  step  is  to  express  At(p)  as  the  sum  of  At^  and  At^  as  depicted 
in  Fig.  II. 8.  Hence, 
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For  large  thickness  ratios  the  travel  time  curves  are  straight  and  pareaiel 
across  dx;  this  implies  that  in  the  thin  layer  case  Atg/Ax  approaches  the 
slope  p.  If  this  substitution  is  made  in  Eq.  11.22,  one  obtains 


At, 


Ax 


+ P 


L ..  ^*1  .2 

"2  - + P 


11.23 


If  Ax  is  eliminated  using  Eq.  11.20,  the  result  is 


2 . 1 

^ U(Az)^ 


II.  2U 


This  equation  does  not  depend  on  Ax  and  At  , which  become  unstable  in 

2 P 

the  case  of  thin  layers.  Equation  II. 2U  is  a linear  equation  in  p and 
2 

At^  and  can  be  solved  graphically  for  and  At.  It  must  be  remembered 


that  At^  is  the  time  difference  between  adjacent  T-X  curves  at  the  same 
value  of  X. 


How  this  method  is  applied  can  best  be  illustrated  by  considering 
an  example.  Figiire  II. 9 depicts  the  travel  time  cvirves  for  a 2-layer 
configuration  in  which  the  first  layer  is  a 5000  m water  layer  with  a 
bilinear  profile  and  the  second  layer  is  a sediment  layer  50  m thick 
having  an  initial  sound  speed  of  1550  m/sec  and  a constant  gradient  of 
0.6  sec  The  bilinear  water  profile  which  was  used  is  listed  on 
Fig.  II. 9 as  a set  of  depths  and  velocities.  The  profile  is  linear 
between  these  points.  The  data  for  the  travel  time  curves  in  Fig.  II. 9 
were  obtained  from  ray  tracing  programs  RANGER  and  PROBE^^  \rtiich  were 
developed  at  ARL/UT  by  Terry  L.  Foreman.  As  is  obvious  from  Fig.  II. 9, 
the  two  curves  are  essentially  parallel  throughout  their  range. 
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The  first  step  in  the  TLA  is  to  calculate  the  ray  parameters 
for  the  top  curve  for  a number  of  values  of  range.  Table  II. 2 shows 
the  X-T  values  of  the  top  curve  which  were  considered  and  the  ray  param- 
eters calculated  ninnerically  with  these  values.  The  next  step  is  to 
measure,  from  the  T-X  cvirves,  the  At^  values  at  the  range  points  considered. 

These  are  also  included  in  Table  II. 2.  The  next  step  is  to  construct 
2 2 

a graph  of  p versus  At^  ; this  graph  is  shown  in  Fig.  II. 10.  The  values 
used  in  the  graph  can  be  fit,  using  least  squares,  to  a straight  line. 

In  the  case  being  considered,  the  least  squares  fit  yielded 

p^  = -102.9**!  At^^  + O.U1682 

where  p is  in  seconds/kilometer  and  At^  in  seconds.  A comparison  with 
Eq.  II. 2U  shows  that 

= O.U1682 

^1 

= 1.5^9  km/sec  , 

eind 

— = 102. 9U1  km~^ 

4(Az)‘^ 

or 


Az  = U9  m 

The  true  value  for  Az  is  50  m and  the  interval  velocity  calctdated  with 
Eq.  A. 6 is  1.565  km/sec.  This  agreement  is  remarkable  considering  the 
thickness  ratio  involved.  But  again  it  must  be  remembered  that  the 
data  being  euialyzed  was  generated  from  a ray  tracing  program.  This 
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TABLE  II. 2 


1 


n 


t 1 


THIN  LAYER  METHOD 


X(tan)  T(sec) 

P(sec/lan) 

At^(sec ) 

1 

6.7096 

0.065 

0.0636 

2 

6.8065 

0.128 

0.0626 

3 

6.9651 

0.188 

0.0610 

1* 

7.1812 

0.2k3 

0.0590 

5 

J.hk99 

0.293 

0.0567 

6 

7.7656 

0.337 

O.05U2 

7 

8.1229 

0.376 

0.0515 

8 

8.5166 

O.UlO 

O.OU89 

9 

8.9»tl8 

O.U39 

O.OU62 

p2  = 

-102. 9i*  + 

0.itl682 

15^9  m/sec 

Az  s 

1*9  m 

li 

■U 
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degree  of  precision  is  almost  never  reeaized  when  using  real  data. 
Bryan  has  determined  that  the  TLA  can  give  good  results  for  thickness 
ratios  of  up  to  500:1. 
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Sound  Speed  Gradient  Information  from  the  RPM* 

In  the  derivation  of  the  RPM  it  is  necessary  to  assume  that 
the  layer  to  be  analyzed  is  Isovelocity.  When  nonhomogeneous  layers 
are  considered,  the  velocity  obtained  by  the  RPM  is  interpreted  as  the 
layer  velocity  emd  no  information  about  the  inhomogeneities  is  obtained. 


There  exists  the  possibility  that  the  RPM  may  be  used  to 

obtain  information  about  sound  speed  gradients  in  a layer  through  the 

use  of  a nomogram  which  gives  possible  combinations  for  C , Az,  and  g. 

o 

The  qujuitity  is  the  velocity  at  the  top  of  the  layer,  Az  is  the  layer 
thickness,  and  g is  the  constant  sound  speed  gradient  in  the  layer,  the 
sound  speed  having  been  assimed  to  be  C(Z)=C^+gz. 


The  RPM  works  with  time  differences  for  rays  traveling  through 
n and  n+1  layers  while  the  range  is  increasing  between  source  and  receiver. 
The  At  between  two  travel  time  curves  changes  with  the  range.  At  range 
zero  (0*90°),  the  At  cam  be  caused  by  an  infinite  number  of  combinations 
of  Az,  c^,  and  g.  There  is,  however,  only  one  combination  of  the  param- 
eters which  fulfills  the  At  values  over  the  entire  range.  In  other  words. 
At  as  a function  of  X is  luiiquely  specified  by  Az,  C^,  and  g. 


The  only  problem  is  that  over  large  ranges  of  X,  travel  time 
curves  belonging  to  a munber  of  different  combinations  of  Az,  C^,  and  g 
are  physicad-ly  Indistinguishable.  This  effect  is  most  noticeable  for 
small  amd  Intermediate  vad.ues  of  range.  The  reason  for  this  is  that  in 
the  relatively  thinner  sediment  layers,  the  amoxmt  of  time  the  ray  spends 
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The  materlaLl  discussed  in  this  subsection  is  a summary  of  previous  research 
done  by  Jan  Drager  as  a visiting  scientist  from  West  Germany. 
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in  the  layer  is  small  ccmipared  to  the  total  time  along  the  ray  path; 
hence  the  travel  time  curve  is  fairly  insensitive  to  the  sediment  param- 
eters. As  the  range  is  increased,  the  time  spent  in  the  sediment  layer 
Increases  and  the  effect  of  the  sediment  psLrameters  becomes  more  noticeable. 
Unfortvmately  at  the  longer  ranges,  the  reflected  arrivals  are  often 
hard  to  differentiate  from  various  other  arrivals  such  as  refracted 
arrivals  ajid  refraction  arrivals  and  mxiltiple  bounce  arrivals. 

As  an  example  of  this  effect,  consider  Fig.  II. 11.  These 
travel  time  curves  were  generated  by  an  isovelocity  water  layer  5 km 
thick  and  a sediment  layer  100  m thick  with  a constant  sound  speed  of 
1280  m/sec.  For  the  remges  considered,  these  curves  are  indistinguishable 
from  those  produced  with  a sediment  layer  having  the  parameters  listed 
in  Table  II. 3,  the  water  leiyer  remaining  the  same.  Of  course  there  are 
a number  of  values  in  between  producing  identical  curves  too. 

The  insensitivity  of  travel  time  curves  to  the  sediment 
parameters  makes  obtaining  information  about  velocity  gradients  a 
difficult  task.  One  approach  to  the  problem  is  to  construct  a family 
of  g versus  curves  for  differing  values  of  layer  thickness.  These 
curves  are  produced  by  ray  tracing  programs  which  use  a water  velocity 
profile  appropriate  to  the  experimental  conditions.  Fig^^re  11.12  is  sm 
example  of  such  a family  of  curves  produced  from  the  ARL  program  PROBE. 

The  water  layer  xised  for  this  example  was  a layer  5000  m thick  having 
a constant  sound  speed  of  1500  m/sec. 

In  principle,  if  one  knew  the  velocity  at  the  top  of  a layer 
and  its  thickness,  one  could  find  g from  a family  of  curves  like  those 
in  Fig.  11.12.  If  the  layer  in  question  is  the  first  sediment  layer, 
this  is  not  unreasonable  since  devices  such  as  ARL's  in  situ  conpressional 
wave  profllometer,^  used  while  taking  shallow  sediment  cores,  measures 
and  the  gradient  in  the  upper  portions  of  the  sediment  layers.  If, 
however,  the  layer  in  question  does  not  sheure  an  interface  with  the  water, 
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TABLE  II. 3 


VALUES  OF  C^,  g,  AND  Az  PRODUCING  IDENTICAL  TRAVEL  TIME  CURVES 


Az(m)  g(sec~^) 

100  lU.5 

5.92 
3.75 
-0- 
-2.8 
-5.0 


C^(at  top  of  layer)  m/sec 

665 

1008 

1100 

1280 

1U26 

1550 
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the  velocity  at  the  top  of  the  layer  is  not  usually  known,  which  in 
turn  makes  it  hard  to  pin  down  the  value  of  g.  It  is  possible,  though, 
in  a case  like  this  to  determine  some  upper  £ind  lower  bound  on  the 
values  of  g.  This  is  done  in  the  following  manner. 

The  thickness  of  a layer  is  usually  known  fairly  well.  In 
the  RPM  it  is  the  parameter  that  is  most  insensitive  to  error.  From 
the  layer  velocity  and  information  about  the  layers  above,  it  is  possible 
to  estimate,  using  physical  intuition,  the  possible  range  of  the  velocity 
at  the  top  of  the  layer.  With  these  bits  of  information,  one  may  consult 
an  appropriate  graph  similar  to  Fig.  11.12  and  determine  a range  of 
possible  values  for  g.  This  information  may  then  be  combined  with  other 
techniques  for  finding  sound  velocity  gradients. 

2 

C.  Comparison  of  the  T (X)  and  Ray  Parameter  Methods 

This  final  section  of  Chapter  II  will  compare  the  two  reflected 
arrival  methods,  the  points  of  comparison  being  accuracy,  ease  of  use, 
and  conditions  under  which  each  method  is  most  applicable. 

When  accuracy  is  considered,  the  RPM  has  the  advantage.  This  is 
because  the  simplifying  assumption  of  isovelocity  layers  is  made  only 
in  the  layer  being  analyzed.  The  velocity  structure  of  the  overlying 
layers  does  not  enter  the  calculations.  In  the  T (X)  method,  the  interval 
velocity  solution  requires  a knowledge  of  all  overlying  layers  and  these 
layers,  which  are  assumed  to  be  isovelocity,  are  added  sources  of  error. 

2 

The  T (X)  method  may  be  improved  when  applied  to  situations  in  which 
the  first  layer  is  a thick  water  layer  with  a known  velocity  profile. 

If  the  water  profile  is  known  then  ray  tracing  techniques  can  be  used 
to  more  accurately  determine  the  amount  of  time  and  range  which  results 
because  of  the  water  layer.  This  added  step  would  greatly  complicate 
the  technique  however. 
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When  one  considers  ease  of  use,  the  T (X)  method  has  the  advantage. 

2 

The  T (X)  method  can  be  used  with  nothing  more  than  a ruler,  graph  paper, 
and  calculator,  whereas  the  RPM  requires  a computer  if  quick  results  are 
needed;  however,  this  is  usually  not  an  important  problem  since  most 
profiling  data  is  not  analyzed  in  the  field. 

2 

The  conditions  under  which  the  T (X)  can  be  used  with  success  are 
those  involving  layers  with  small  thickness  ratios.  Bryan^  has  determined 
that  for  thickness  ratios  less  than  15  the  interval  velocity  in  the  thin 
layer  is  not  sensitive  to  the  velocity  structtire  above.  If  one  is  dealing 
with  small  thickness  ratios  (shallow  water,  for  example),  then  the  T (X) 
methods  give  good  results  quickly.  Also,  the  T (X)  method  as  developed 
by  Durbaiam  is  useful  when  sloping  subbottom  layers  are  considered. 


The  RPM  can  be  applied  to  a wide  range  of  layer  conf igtirations . 

The  normal  RPM  can  be  used  for  thickness  ratios  up  to  50:1.  When 
thinner  layers  are  encountered,  the  thin  leiyer  approximation  may  be 
used  out  to  thickness  ratios  of  500:1.  The  ability  to  handle  thin 
layers  makes  the  RPM  most  useful  when  analyzing  high  resolution  profiling 
data  in  which  many  thin  reflecting  horizons  are  observed. 

Another  condition  for  the  RPM  is  that  relatively  continuous  travel 
time  data  be  available.  This  type  of  data  is  needed  because  a relatively 
fine  mesh  of  T and  X values  is  needed  in  order  to  calculate  the  deriva- 
tives required  by  the  RPM. 


III.  ANALYSIS  OF  INTERVAL  VELOCITY  DATA 


In  the  study  of  sound  propagation  in  the  ocean,  whether  ray  theory 
or  normal  mode  theory  is  used,  one  needs  to  detail,  as  closely  as  ; 

possible,  the  subbottom  velocity  profile.  Unfortunately,  the  subbottom  i 
is  the  region  which  is  the  hardest  to  obtain  information  about.  With  i 

I ; 

the  exception  of  coring  data,  the  layer  and  velocity  structure  of  the  i 

subbottom  must  come  from  the  analysis  of  profiling  data. 

T 0 

The  method  developed  by  Houtz  et  al.  ’ will  allow  one  to  obtain 
the  velocity  profile  of  the  subbottom  from  interval  velocities  calculated 
using  the  methods  of  Chapter  II.  For  lack  of  a better  name,  this  method  3 
will  be  referred  to  as  the  V(t)  method,  since  it  works  with  velocity  as 
a function  of  one-way  travel  time.  The  small  letter  "t"  will  be  used  ^ 

in  the  chapter  to  denote  one-way  travel  time  as  opposed  to  a capital  "T" 
used  in  the  previous  chapter  to  denote  two-way  time.  1 

A.  Derivation  of  the  V(t)  Method 

The  V(t)  method  is  one  in  which  interval  velocities  and  travel  ! 

times  axe  used  to  construct  an  instantaneous  velocity  fxuiction,  V(t), 
from  which  a gradient  and  a depth  function  can  be  obtained.  The  first 
step  in  the  derivation  of  this  method  is  to  relate  the  interval  velocities 
to  the  instantaneous  velocities.  In  the  V(t)  method,  the  interval  velocity 
of  a layer  is  assumed  to  be  the  instantaneous  velocity  at  the  midpoint, 
in  time,  of  that  layer.  Or,  in  other  words,  if  a certain  reflecting  ^ 

layer  is  bounded  by  one-way  times  t and  t+At,  then  the  interval  velocity 
of  the  layer  is  assmed  to  be  the  instantaneous  velocity  at  t+At/2.  The 
assignment  of  the  interval  velocity  to  the  instantaneous  velocity  at 
the  midpoint  in  time  is  not  entirely  correct  but  Hc^tz  et  al.  in  their  i 

paper,  "Velocity  of  Deep-Sea  Sediments  from  Sonobouy  Data,"  show  that  I 

this  assumption  is  generally  in  error  by  much  less  than  50  m/sec.  j 


The  V(t)  meohod  is  implemented  by  fitting  the  interval  velocities 
to  a continuous  third  order  polynomial  in  t,  the  one-way  travel  time,  i.e.. 


V(t)  = A + Bt  + Ct^  + Dt^ 


III.l 


A continuous  velocity  function  is  used  in  Eq.  III.l,  even  though  velocity 
discontinuities  are  quite  common,  because  one  is  dealing  with  instantaneous 
velocities  at  the  midpoints  of  layers.  Profiling  techniques  do  not  provide 
enough  information  to  handle  velocity  discontinuities  which  would  occur 
at  the  layer  boundaries  between  adjacent  interval  velocities;  the  effect 
is  that  any  velocity  discontinuities  that  might  be  present  are  smoothed 
out.  An  example  of  this  effect  will  be  shown  later  in  this  chapter. 


A depth  function  may  be  obtained  from  Eq.  III.l  by  integrating  with 
respect  to  t to  give 


H(t)  = At  + I Bt^  + j Ct^  + ^ Dt** 


III.  2 


A gradient  function  may  be  obtained  from  Eq.  III.l  in  the  following 
manner.  If  V(t)  is  written  as 


v(0  = ^ 


then  the  gradient  as  a function  of  t is  formed  from  the  following  ratios 


g(t)  = 


<iV(t) 

dt 

dH(t) 

dt 


dv(t) 

dt 

v{t) 


, . X B + 2Ct  + 3Dt 

git;  = = ^ 

A + Bt  + Ct  + Dt^ 


Once  the  coefficients  A,  B,  C,  and  D sure  known,  the  depth  and  gradient 
as  a function  of  t are  easy  to  calculate. 

If  a layer  is  bounded  by  one-way  times  t and  t+At,  then  the  layer 
thickness,  AZ,  may  be  found  using  the  following, 

AZ  = H(t+At)  - H(t)  . ] 

The  gradient  of  this  same  layer  may  be  reasonably  well  approximated  by 
the  gradient  function  evaluated  at  the  midpoint  in  time  of  this  layer. 


g - g(t  + ^ At)  . 


If  a continuous  velocity  profile  through  the  subbottom  is  desired,  i.e., 
V as  a function  of  H,  it  may  be  obtained  from  Eqs.  III.l  suid  III. 2 by 
forming  pairs  of  velocity  smd  depth  vsilues  for  the  same  values  of  t. 


At  ARL/UT,  a slightly  different  V(t)  method  has  been  used.  A 
computer  program  written  by  Terry  Foremsui  fits  the  interval  velocities 
and  the  times  using  the  cubic  spline  technique, rather  than  using 
least-squares  to  fit  them  to  a third  order  polynomial.  The  cubic  spline 
method  is  one  in  which  a cubic  equation  like  Eq.  III.l  is  fit  between 
each  pair  of  velocity-time  values.  Instead  of  one  set  of  fitting 
coefficients,  there  are  a nxunber  of  sets  of  fitting  coefficients,  each 
set  being  used  between  a particular  pair  of  adjacent  velocity-time 
values.  To  evaluate  the  sets  of  coefficients,  each  cubic  spline  is 
required  to  match  the  adjacent  splines  with  continuous  first  and  second 
derivatives.  The  end  result  of  this  procedure  is  a velocity  function 
which  is  smoother  eind  has  less  abrupt  behavior  than  one  obtained  by 
least  squares. 


B. 

As  an  example  of  the  V(t)  method,  consider  the  case  of  four  subbottom 
layers  xinder lying  a 5000  m water  layer  having  a bilinear  profile.  Table  III.l 
lists  the  parameters  of  the  subbottom  layers  and  the  water  layer.  The 
subbottom  profile  of  this  example  has  gradient  discontinuities  across  the 
layer  interfaces  but  the  sound  velocities  are  continuous. 

Figure  III.l  depicts  the  computer  generated  T-X  curves  for  the 
example.  Since  the  thickness  ratios  are  large,  the  ray  parameter  method, 
using  the  thin  layer  approximation,  was  employed  to  determine  the  interval 
velocities,  layer  thicknesses,  and  vertical  travel  times.  These  values 
are  listed  in  Table  III. 2.  For  comparison  sake,  the  actual  layer  thick- 
nesses and  the  layer  velocities  calculated  from  Eq.  A. 6 are  included  in 
parenthesis. 

Before  applying  the  V(t)  method,  one  needs  to  form  a set  of  values 
for  the  Interval  velocity  versus  one-way  travel  time  to  the  midpoints, 
in  time,  of  the  four  subbottom  layers.  These  values  are  listed  in 
Table  III. 3.  The  one-way  times  are  referenced  to  the  interface  between 
the  water  and  the  first  subbottom  layer,  i.e.,  t*0  refers  to  the  top  of 
the  first  subbottom  layer. 

The  set  of  velocity-time  values  is  next  used  as  input  to  the  ARL/UT 
spline  fit  program.  This  program  determines  the  instanteuieous  velocity 
function,  the  depth  function,  and  the  gradient  function  euid  prints  them 
in  table  form  for  a fine  mesh  of  t values.  In  this  example,  the  value 
of  instantaneous  velocity  at  the  top  of  the  first  layer,  i.e.,  t*0,  is 
assumed  to  be  known  and  is  included  in  the  fitting  procedure.  The 
program  also  generates  Calcomp  plots  of  the  instantaneous  velocity 
function  and  the  velocity  as  a function  of  depth.  These  plots  are 
given  in  Figs.  III. 2 and  III. 3. 
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TABLE  III.l 


Water  Profile 


Depth 

(m) 


Sound  Velocity 
(m/sec ) 


Subbottom  Profile 


Subbottom 
Layer  No. 

Thickness 

(m) 

Constant 

Gradient 

{sec”^) 

Initial  Sound 
Velocity 
(m/sec ) 

1 

100 

0.5 

1550 

2 

100 

1.0 

1600 

3 

100 

1.5 

1700 

1* 

100 

2.0 

1850 

TABLE  III. 2 


r 


LAYER  PARAMETERS  CALCULATED  BY  THE  THIN  LAYER 
APPROXIMATION  OF  THE  RAY  PARAMETER  METHOD 


Subbottom 
Layer  No. 

AZ 

lai 

Interval  Velocity 
(m/sec ) 

2-Way  Vertical 
Travel  Time 
(sec) 

1 

99  (100) 

1560  (1575) 

0.1270 

2 

101  (100) 

1660  (1650) 

0.1212 

3 

100  (100) 

' 1779  (I77U) 

0.1127 

U 

101  (100) 

1957  il9h9) 

0.1027 

TABLE  III. 3 

INTERVAL  VELOCITY  AS  A FUNCTION  OF  ONE-WAY  TIME 

One-Way  Time  to 

Interval  Velocity  Midpoint  of  Layer 

(m/sec) (sec) 

0.03175 
0.09381 
0.15228 
0.20613 
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1560 
1660 
1779 

L'1  1957 


WATER  PROFILE 


X - km 


J 

12 


FIGURE  m.1 

TRAVEL  TIME  CURVES  FOR  A THICK  WATER  LAYER 
OVERLYING  FOUR  THIN  SEDIMENT  LAYERS 
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FIGURE  m.2 

INSTANTANEOUS  VELOCITY  AS  A FUNCTION  OF  ONE-WAY 
TIME  OBTAINED  FROM  INTERVAL  VELOCITIES 


ARL  . UT 
AS- 76 -1113 
SRR  • DR 
11-2-76 


From  the  table  of  velocity,  depth,  and  gradient  values  it  is 
possible  to  reconstruct  the  subbottom  profile.  The  layer  thickness  and 
layer  gradient  may  be  approximated  by  Eqs.  IIl.U  and  III. 5,  while  the 
velocity  at  the  top  of  each  layer  may  be  approximated  by  the  instantsmeous 
velocity  function  evaluated  at  the  one-way  time  to  the  top  of  the  layer 
in  question.  The  velocity  profile  determined  in  this  manner  for  the 
present  example  is  listed  in  Table  III.U.  Again,  the  correct  values 
are  listed  in  parenthesis. 

Table  III.U  shows  that  the  V(t)  method  gives  very  good  results 
whenever  good  interval  velocity  values  ceui  be  obtained  and  when  there 
are  no  velocity  discontinuities  between  the  layers.  Unf ortionately , 
velocity  discontinuities  are  common  and  good  interval  velocity  values 
are  frequently  hard  to  obtain  when  analyzing  real  data.  These  points 
will  be  addressed  further  in  the  next  section. 

A comparison  of  Tables  III. 2 and  III.U  shows  that  the  values  of 
layer  thickness  obtained  by  the  ray  parameter  method  seem  to  be  better 
than  those  obtained  by  the  V(t)  method.  The  RPM  generally  yields  a good 
value  for  layer  thickness.  By  combining  these  two  methods,  even  a better 
approximation  to  the  real  profile  may  be  obtained. 

C.  Problems  Associated  with  the  V(t)  Method 

This  section  discusses  problems  which  arise  when  the  V(t)  method 
is  being  employed,  the  problem  due  to  velocity  discontinuities  and  the 
effect  of  measurement  errors  on  the  results  of  the  V(t)  method. 

1.  Velocity  Discontinuities 

As  seen  in  the  example  of  Section  B,  when  no  velocity 
discontinuities  are  present  the  V(t)  method  is  capable  of  reconstructing 
the  subbottom  velocity  profile  quite  well.  When  velocity  discontinuities 
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TABLE  III.U 


, 1 


i ' 


VELOCITY  PROFILE  RECONSTRUCTED  FROM  V(t)  METHOD 
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Subbottom 
Layer  No. 

Thickness 

(m) 

Consteuit 

Gradient 

(sec“^) 

Initial  Sound 
Velocity 
(m/sec) 

1 

99.7  (100) 

0.5  (0.5) 

1550  (1550) 

2 

95.2  (100) 

1.1  (1.0) 

1602  (l600) 

3 

77.0  (100) 

1.5  (1.5) 

1716  (1700) 

k 

129.0  (100) 

1.8  (2.0) 

1862  (1850) 

are  present,  certain  problems  in  the  V(t)  method  arise.  These  problems 
are  best  illustrated  by  considering  two  examples. 


Consider  two  configurations,  each  having  four  subbottom  layers. 

The  velocity  profiles  of  these  subbottom  layers  are  listed  in  Table  III. 5. 

The  first  configuration  is  one  in  which  there  are  no  velocity  discontinuities 
between  the  layers.  The  second  configuration  is  identical  to  the  first 
with  the  exception  that  velocity  discontinuities  between  the  layers  are 
allowed. 

For  the  purposes  of  this  illustration,  it  was  not  necessary  to 
analyze  travel  time  curves.  Instead,  the  interval  velocities  and  travel 
times  to  be  used  in  the  V(t)  method  were  calculated  using  the  equations 
of  Appendix  A.  These  calculated  veulues  are  shown  in  Table  III. 6.  As  in 
the  last  example,  these  values  were  used  in  the  spline  fit  program  to 
compute  instantaneous  velocity,  depth,  and  gradient  functions.  The  sub- 
bottom profiles  reconstructed  from  the  results  of  this  program  are  listed 
in  Table  III. 7. 

A comparison  of  Tables  III. 5 and  III. 7 shows  that,  in  the 
continuous  case,  the  V(t)  method  reproduces  the  subbottom  profile  quite 
well,  with  the  exception  of  the  thicknesses  of  the  last  two  layers  which 
are  in  error  by  about  13^.  The  gradient  values  and  the  Interval  sound 
velocities  were  reproduced  quite  accurately. 

The  discontinuous  case,  however,  did  not  work  as  well.  The 
errors  in  the  thickness  values  and  in  the  Initial  soxind  speed  all 
increased  in  the  discontinuous  case,  thotigh  the  initial  sound  speed 
proved  to  be  less  sensitive  (maximum  difference  of  1.5JS).  The  most 
pronounced  differences  were  observed  in  the  predicted  gradients,  the 
average  percent  difference  being  30%.  The  reason  for  this  difference 
is  easily  explained.  As  mentioned  before,  the  V(t)  method  obtains  a 
continuous  instantaneous  velocity  function.  Whenever  velocity 


58 


I 

i 


TABLE  III. 5 


Subbottom 
Layer  No. 


1 

2 

3 

U 

1 

2 

3 

U 


SUBBOTTOM  PROFILES 

FOR  EXAMPLES  OF  SECTION 

III.C 

Thickness 

(m) 

Constant 

Gradient 

Initial  Sound 
Speed 

(sec“^) 

(m/sec) 

Continuous  Case 

100 

1 

1550 

100 

1 

1650 

100 

1 

1750 

100 

1 

1850 

Discontinuous  Case 

100 

1 

1550 

100 

1 

1700 

100 

1 

1820 

100 

1 

i960 
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TABLE  III. 6 

I 


LAYER  PARAMETERS  CALCULATED  USING  THE  RESULTS  OF  APPENDIX  A 


Subbottom 
Layer  No, 

Two-Way  Vertical 
Travel  Time 
(sec) 

Interval 
Velocity 
(m/sec ) 

One-Way 
Midpoint  Time 
(sec) 

Continuous  Case 

1 

0.1250 

1600 

0.03126 

2 

0.1177 

1700 

0.09192 

3 

0.1111 

1800 

O.IU913 

k 

0.1053 

1900 

0.20322 

Discontinuous  Case 

1 

0.1250 

1600 

0.03126 

2 

O.IIU3 

1750 

0.09108 

3 

0.1070 

1870 

O.IU639 

U 

0.0995 

2009 

0.19803 

TABLE  III. 7 

VELOCITY  PROFILES  RECONSTRUCTED  FROM  THEORETICAL 
VALUES  OF  INTERVAL  VELOCITIES  AND  TRAVEL  TIMES 


Continuous  Case 


1 

98.5 

1.006 

1550 

2 

95.5 

0.996 

1651 

3 

87.1 

l.OOU 

1751 

1* 

87.1 

0.98k 

1851 

Discontinuous  Case 


1 

98. U 

0.905 

1550 

2 

93. U 

1.379 

1675 

3 

73.2 

I.2U9 

1812 

U 

119.  U 

I.U29 

1939 

SubUottom  Thickness 

Layer  No.  (m) 


Constant 

Gradient 

(sec~^) 


Initial  Sound 
Speed 
(m/sec) 


discontinuities  are  present,  the  interval  velocity  of  a layer  is  larger 
than  it  would  be  if  no  discontinuities  were  present.  Since  the  V(t) 
method  does  not  recognize  discontinuities,  it  reacts  by  predicting  a 
larger  value  for  the  layer  gradient  to  compensate  for  the  higher 
interval  velocity. 

This  brings  up  the  question  of  whether  velocity  discontinuities 
can  be  detected  using  the  V(t)  method.  The  answer  has  to  be,  "sometimes.” 
Anytime  the  V(t)  method  gives  gradient  values  that  seem  vuirealistically 
high,  the  expected  range  being  0.5  sec“-^  to  2.0  sec"-^,  one  can  suspect 
that  a velocity  discontinuity  exists.  In  general,  to  determine  the 
gradient  more  closely,  one  must  use  and  combine  as  much  infonnation  and 
as  many  methods  as  possible.  For  example,  if  velocity  discontinuities 
are  suspected,  the  V(t)  method  can  be  used  to  produce  a 'first  guess.' 

Well  cores  from  geologically  similar  regions  and  possibly  the  gradient 
analysis  discussed  in  Chapter  II  can  then  be  used  to  revise  the  first 
guess.  In  this  manner  one  can  obtain  a best  possible  guess  of  the 
velocity  profile. 

2.  Effects  of  Measurement  Errors 

All  of  the  data  of  the  examples  considered  in  this  chapter 
were  generated  from  a computer  ray  tracing  program.  It  goes  without 
saying  that  Interval  velocities  calculated  from  actual  data  will  have 
significant  uncertainties  associated  with  them,  and  these  \uicertainties 
will  affect  the  subbottom  velocity  profiles. 

13  2 

Clay  and  Rona  indicate  that  the  T (X)  method  of  intervad. 
velocity  analysis  produces  \incertainties  of  the  order  of  ±200  m/sec. 

It  is  expected  that  this  uncertainty  can  be  reduced  to  the  order  of 
±100  m/sec  by  use  of  the  ray  parameter  methods. 
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To  get  a feeling  for  the  effect  this  uncertainty  has  on  the 
gradient,  consider  the  following  simplified  'linear  analysis.'  Assume 
that  one  is  attempting  to  determine  the  gradient  in  a single  sediment 


layer  having  an  initial  sound  speed  of  1500  m/sec,  a thickness  of  300  m, 
and  a constant  gradient  of  1 sec  The  equations  derived  in  Appendix  A 


yield  an  interval  velocity  value  of  l6U8  m/sec  and  a two-way  vertical 
travel  time  of  0.36U6  sec,  or  a one-way  time  to  the  midpoint  of  0.0912  sec 
(see  Fig.  III.U). 


The  depth  coordinate  of  this  interval  velocity  may  be 
approximated  by  the  product  of  intei*val  velocity  and  one-way  midpoint 
time  [Z=(l6W  in. /sec) (0.0912  sec)  = 150.3  m]. 


The  simplified  linear  analysis  mentioned  above  consists  of 
calculating  the  linear  slope  between  the  two  velocity  depth  values, 
(1500  m/sec,  0)  and  {16U8  m/sec,  150.3  m).  The  result  is 


=,  . 0.98  seo-^ 

m 


If  the  1500  m/sec  velocity  is  assumed  to  be  exact,  suid  uncertainties 
of  ±200  m/sec  and  ±5  m are  assumed  for  the  interval  velocity  and  the 
depth  coordinate,  statistical  analysis  for  the  combination  of  errors 


for  the  above  eqiiation  yields  an  vincertainty  for  g of  about  ±1.1*  sec 


-1 


If  the  vincertainty  of  the  Interval  velocity  is  taken  to  be  ±100  msec, 

^ Hence,  errors  in  the 


the  uncertainty  in  g is  reduced  to  ±0.8  sec 
measurement  of  interval  velocities  have  a tremendous  effect  on  the 
confidence  one  has  in  his  results. 


To  see  the  possible  effects  of  measurement  errors  on  a 
multilayer  case,  consider  Fig.  III. 3.  In  this  figure,  the  conservative 
value  for  the  uncertainty  of  the  interval  velocities,  i.e.,  ±100  m/sec, 
is  marked  by  the  vertical  bars  centered  on  the  interval  velocity  points 
calculated  by  the  RPM.  For  comparison  the  Interval  velocities  calculated 


IV.  REFRACTED  ARRIVAL  METHODS 


The  two  methods  to  be  discussed  in  this  chapter  obtain  the  velocity 
profile  in  a single  subbottom  layer  by  techniques  that  employ  the  effects 

Q 

of  refracted  arrivals.  The  first  method,  used  by  Dicus,  works  with 
travel  times  of  the  refracted  arrivals.  The  second  method,  developed  by 
Hanna, calculates  transmission  loss  due  to  energy  refracted  through 
the  subbottom  as  a function  of  range. 


Both  these  methods  use  the  same  model  for  the  subbottom.  The  assvimed 
bottom  configuration  is  a single  layer  which  has  a pseudolinear  velocity 
profile  given  by 


where  is  the  velocity  at  the  top  of  the  layer  eind  3 is  the  sound  speed 
gradient  at  the  surface  of  the  layer  (3=C'(0)). 


A.  Dicus'  Method 


Dicus'  method  is  capable  of  determining  the  initial  values  of  sound 
speed  and  gradient  in  a single  layer  through  an  anaQ.ysis  of  the  time 
differences  berween  arrivals  that  refract  through  the  bottom  and  those 
that  reflect  from  the  top  of  the  subbottom  layer  (see  Fig.  IV.l).  Dicus' 
method  is  derived  in  the  following  manner. 


1.  Derivation 


From  Fig.  IV.l,  one  may  write  the  following  expression  for  the 
time  difference  between  paths  1 and  2, 


C(z)  I C(^ 
path  1 •'path  2 

The  quantity  ds  is  an  incremental  distance  along  the  ray  path.  If  the 
water  layer  is  assumed  to  have  a slowly  varying  velocity  profile,  then 
the  water  paths  may  be  approximated  by  straight  lines  and  Eq.  IV. 2 becomes 

/2S 

’ IV.3 

i-cith  1 

where  is  the  time  averaged  sound  speed  in  the  water. 

From  the  geometry  of  Fig.  IV. 1,  the  integral  in  Eq.  IV. 3 may 
be  written  as 


L 


path  1 


ds 

cM 


IV. 


where  sed.  refers  to  that  portion  of  path  1 which  lies  in  the  sediment 
layer.  Ccmbining  Eqs.  IV. 3 and  IV. U,  one  obtains 


At 


ds  2AS 


IV. 5 


with  ASsSg-S^.  If  AS  is  approximated  by  AX  cos  (|),  Eq.  IV. 5 becomes 


ds  2AX  cos  d) 


The  integral  In  Eq.  IV. 6 may  be  written  as 


/ ctIt’ 

•'sed.  •'o 


AZ 


dZ 

sin  e(z)  C(Z) 


IV. 7 


where  dZ=ds  sin  0(Z)  has  been  used.  The  second  integral  in  Eq.  IV. 7 is 
most  easily  evaluated  by  changing  the  variable  of  Integration  from  Z to 
C{Z),  the  sound  velocity  in  the  sediment  layer.  This  procedure  yields 


/ 

•sed. 


, 2C 

ds o 

C(Z)  6 cos  0 


/ 


C /cos  0 
o o 


dC 


C^(Z) 


° - C^(Z) 


cos^  0 


' 1 


' I 

f 1 


2 sin  0 
c 

30 


^1+2  cos^  0^^ 


To  obtain  Eq.  IV. 8,  the  following  relationships  were  used, 

8 


dZ 


Bc^(z) 


sin  0(Z)»  Jl  - cos^  0(Z)  * /l  - C‘"{Z) 


» cos^  0 
\ o 


IV.  8 


If  Eq.  IV. 8 is  substituted  into  Eq.  IV. 6 the  following  expression  is 
obtained , 

2 sin  0 , 


Ax  = 


30 


e \ . 2_c^ 


1+2  cos  0 
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AX 


IV. 9 


To  finish  the  evaluation  of  At,  an  expression  for  AX  in  terms 
of  the  sediment  parameters  is  required.  This  expression  is  obtained  as 
follows , 


/AZ  -C  /cos  0 _ 3^_ 

dZ  _ / ° ° ^o 

tan  0(Z)  ~ I \ 3,  ^ 

tan  0(Z)0  C-^(Z) 


If  one  uses 


C^(Z)  cos^  0 


tan  0(Z)  = 

cos  9(Z) 


C(Z)  . 
— = — cos  0 

^o 


in  Eq.  IV. 10,  the  result  is 


3 rC  /cos  0 
/ o o 


■¥/ 


C cos  0 sin  0 
o o o 


c y-i—  - 

Vcos  0 

o 


With  Eq.  IV. 11,  Eq.  IV. 9 becomes 


2 sin  0 / ? \ o A 0 sin  0 

AX  =-3^(i  + 2cos  gO o 


If  the  water  layer  is  thick  and  the  source  depth  shallow,  the  euigle  4> 
may  be  approximated  well  by  the  angle  0w.  This  substitution  in 
Eq.  IV. 12  yields 
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T 


2 sin  6 


At  =s 


, „ \ 2 cos  6 C cos  0 sin  0 

f a , 


From  Snell's  Law,  cos  ®ay  ^ replaced  by  cos  0^  to  give 


2 sin  0 


At  * 


r.  f 9 \ 9 

(i  1 2 =0=2  ej  - 2 COB®  — g-a 


(l  - COB®  6^)  = ^ (i  - 


2 sin  0 . p 

— sr^  1-=°=  ®o 


a 0 

cos  0 I 


IV.  Ik 


Snell's  Law  may  be  used  once  again  to  yield 

.3/2 


2 f ^ ''o  2ft 

^w 


IV.  15 


Equation  IV. 15  may  be  transformed  into  a form  that  can  be 
solved  graphically  for  Cq/C^  and  0 by  making  the  following  substitutions. 
Let  r*cos^  0 and  s=(3At/2)2/3,  With  these  changes  of  variables. 


Eq.  IV. 15  can  be  written  as 


,2/3 


1/c 

fe) ' * ^ 

• 

IV.  16 

1 \ w/ 

- 

Equation  IV. l6  Is  a linear  equation  in  s and  r which  can  be  solved  for 
0 and  C^/C^  graphically. 

It  is  possible  to  obtain  the  simple  expression  in  Eq.  IV. l6 
because  of  the  assimption  of  a pseudolinear  profile.  The  pseudolinear 
profile  does,  however,  have  a few  problems  associated  with  it.  At  a 
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critical  depth  of  pseudolinear  profile  is  singiaar  and 

as  Z approaches  'the  velocity  profile  becomes  physically  unrealistic 

This  can  have  an  effect  on  Dicus’  method  if  one  is  analyzing  data  from 
thick  subbottom  layers.  For  this  reason,  when  thick  layers  are  present 
one  should  analyze  that  part  of  the  data  which  has  the  smaller  bottom 
angles,  i.e.,  the  shallower  refracting  rays. 


If  thin  layers  or  shallow  refracted  rays  are  considered,  the 
pseudolinear  gradient  parameter  6 should  be  nearly  the  same  as  the 
constant  gradient  g.  From  Eq.  IV. 1 in  the  limit  of  small  values  of  the 
ratio  2BZ/C  , one  may  write 


Hence,  «.n  tula  limit,  H should  be  very  nearly  equal  to  g,  the  linear 
gr^.41e(it 


«tf;.ie  of  Dicua*  Method 


io  llTvet  »*•  this  method,  consider  a 2-layer  configuration 
like  the  ct  * .n  Fig.  IV.  1.  Let  the  top  layer  correspond  to  an 

I ipy  water  lajrer  with  a sound  speed  (C^)  of  1500  n/sec  and 
thick...  :i  -f  5000  a.  Let  the  subbotton  layer  be  a layer  UOO  m thick 
wit  .*.<  initial  souKf  velocity  of  1550  m/sec  and  a constant  gradient 
of  g*l.C  «ec  (The  linear  gradient  is  used  in  this  example  to  show 
the  rorreepondence  between  0 and  g and  because  the  ray  tracing  programs 
asT>(t  to  generate  the  data  of  this  example  require  that  constant  gradients 
be  used. ) 


Table  IV. 1 lists  the  values  of  6^,  At,  s,  and  r which  were 
obtained  from  the  ARL  ray  tracing  program  RANGER.  To  simulate  actual 


'BLE  IV. 1 


QUAMTITIES  USED  IN  THE  APPLICATION  OF  DICUS'  METHOD 


Range 

e 

w 

(deg) 

Ax 

(sec ) 

s . (1  4x) 

(sect's) 

2 

r ■ cos 

15 

37.999 

0.1310 

0.3380 

0.6210 

16 

35.1+98 

0.1051 

0.2918 

0.6628 

17 

33.331+ 

0.0850 

0.2531 

0.6980 

18 

31.1+31+ 

0.0690 

0.2205 

0.7280 

19 

29.71+9 

0.0563 

0.1925 

0.7538 

20 

28.2I1O 

0. u461 

0.1685 

0.7761 

21 

26.879 

0.0378 

0.1476 

0.7956 

22 

25.61+6 

0.0311 

0.1296 

0.8127 

23 

24.521 

0.0255 

0.1135 

0.8277 

2U 

23.492 

0.0209 

0.0994 

0.8411 

s » -1.0833  r + 1.0097‘» 
s = 1.0097U(-1.0728Ur+l) 

= 1.0358 

w 

8 = 0.99  sec”^ 


■1 

! 


experimental  conditions,  the  angle  was  used  to  approximate  0 . 

Figure  IV. 2 depicts  the  graph  of  s versus  r.  From  the  very  nearly 
linear  relation  between  s and  r,  it  would  appear  that  in  this  case  the 
pseudolinear  profile  well  approximates  the  linear  profile  and  that 
replacing  0^  by  ()»  is  Justified. 

A least-squares  linear  fit  between  s and  r yields  the  following 

equation, 

s = 1.0097i*(-1.0T28J+r+l) 

The  values  of  C^/C^  and  8 may  be  obtained  from  the  above  equation  as 
follows , 


7^  = ^0728U  = 1.0358 


w 


3 = 


(1.0097*+)^^^ 


= 0.99  sec 


-1 


The  actual  values  for  these  parameters  are 


“ 1550/1500  = 1.0333  and  g = 1 sec”^ 

The  agreement  between  3 and  g indicates  that  the  linear  profile  used  to 
generate  the  data  is  well  approximated  by  a pseudolinear  profile. 

3.  Discussion 


To  be  able  to  apply  this  method  to  experimental  data  two 
conditions  must  hold.  First,  it  is  necessary  that  one  be  able  to 
identify  refracted  arrivals  in  the  data.  Sometimes  this  is  not  an 
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easy  task.  Usually,  other  preprocessing  such  as  deconvolution  processing 
is  required  before  refracted  arrivals  can  be  identified.  The  second 
condition  for  the  success  of  thj.s  model  is  that  the  angle  0 be  very 
nearly  equal  to  ())  (see  Fig.  IV. 1).  This  is  necessary  because  the  angle  9 
cannot  be  determined  beforehand  without  a knowledge  of  the  subbottom 
profile.  For  this  reason  it  is  necessary  to  approximate  9 with  the 

V 

angle  (j),  which  can  be  determined  if  the  water  profile  and  horizontal 
range  are  known.  In  situations  involving  thick  water  layers  and  shallow 
sources  this  approximation  is  usually  Justified. 

B.  Hanna's  Method 


Hanna's  Method  is  a \inique  method  for  determining  subbottom  velocity 
structiires  by  analyzing  nulls  in  a transmission  loss  curve  for  a bottomed 
receiver.  This  technique  is  described  in  J.  S.  Hanna's  paper,  "Short-Range 
Transmission  Loss  and  the  Evidence  for  Bottom-Refracted  Energy. This 
section  is  essentially  a brief  summary  of  Hanna's  paper. 

1.  Proposed  Method  of  Solution 


'-f ; 


Hanna's  technique  proposes  to  determine  sediment  velocity  and 
gradient  infonnation  from  the  interference  structure  exhibited  by  the 
transmission  loss  curve  for  a bottomed  receiver  over  the  range  of  the 
direct  arrival.  The  range  is  restricted  in  this  manner  because,  outside 
the  range  of  the  direct  arrival,  multibounce  ray  paths  further  complicate 
the  transmission  loss  by  adding  more  interference  structure. 


If  there  were  no  energy  being  returned  from  the  bottom,  the 
.nterference  structure  of  the  transmission  loss  should  be  accounted  for 
1«it*  well  by  the  Uoyd's  Mirror  Effect.  (The  Lloyd's  Mirror  Effect  is 
aotlcMblc  when  a shallow  source  is  used  and  results  from  the  inter- 


ray paths  that  are  upgoing  and  downgoing  at  the  source.) 

allowed  to  refract  or  reflect  back  into  the  water  column. 


1 

^1 
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>: 
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the  transmission  loss  curve  will  show  added  interference  structure. 

Hanna's  method  proposes  to  identify  interference  nulls  in  transmission 
loss  curves  that  arise  because  of  refracted  energy  and  to  determine 
velocity  and  gradient  information  from  the  positions  of  these  nulls. 

2.  Discussion  of  the  Technique 

The  first  assumption  made  when  applying  this  method  is  that 
the  energy  incident  on  the  bottom  can  be  approximated  by  a plane  wave 
and  that  a plane  wave  reflection  coefficient  may  be  used  to  describe 
the  subbottom.  Let  the  reflection  coefficient  be  denoted  by  R(0). 

This  plane  wave  reflection  coefficient  is  specified  as  a function  of 
grazing  angle  and  includes  the  effects  of  energy  refracted  in  the  bottom. 

If  a bottomed  hydrophone  is  used  to  measure  pressure  levels, 
the  amplitude  as  a function  of  grazing  angle,  a(0),  may  be  specified  as 

a(0)  = A(0)[1+R(0)]  . IV. 17 

The  function  A(0)  is  the  amplitude  of  the  incident  plane  wave  and  contains 
information  concerning  spreading  loss  and  Lloyd's  Mirror  Effect.  The 
term  1+R(0)  contains  information  about  interference  between  incident  and 
reflected  or  refracted  energy. 

Before  this  method  can  be  developed  further,  a model  for  the 
water  and  sediment  layers  must  be  specified.  Hanna  models  the  water 
layer  as  isovelocity  and  the  subbottom  as  a sediment  layer  with  a 
constant  gradient  overlying  a solid  isovelocity  basement  (see  Fig.  IV. 3). 
For  this  configuration,  the  plane  wave  reflection  coefficient  may  be 
approximated  by  the  reflection  coefficient  for  plane  waves  incident  from 
an  isovelocity  haJ-f-space  on  the  boundewy  of  another  half-space.  This 
approximation  should  be  good  for  thick  sediment  layers.  The  reflection 
coefficient  for  this  case  may  be  written  as 
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FIGURE  IZ.3 

LAYER  CONFIGURATION  AND  VELOCITY  STRUCTURE 
(AFTER  HANNA) 
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ik  sin  0iMo)  - (p/p^^)  i^'(O) 


i 


IV.  18 


The  variable  k is  the  wave  number  in  the  isovelocity  top  (water)  layer 
(k=2TTf/c,  c=velocity  in  top  layer),  p and  p^  are  the  densities  in  the 
upper  and  lower  half-spaces,  and  iJi(O)  is  the  wave  function  of  the  lower 
half-space  evaluated  at  the  boundary  between  the  upper,  isovelocity 
half -space. 

The  nulls  in  transmission  loss  curves  which  are  caused  by 
refracted  energy  are  determined  by  the  zeros  of  1+R(0).  From  Eq.  IV. l8, 
this  condition  implies 


sin  0=0 


iKO)  = 0 


IV.  19 


IV.  20 


Equation  IV. 19  is  satisfied  when  the  direct  path  attains  grazing  incidence 
on  the  bottom  (0=0).  The  nulls  specified  by  Eq.  IV. 20  depend  on  the  wave 
function  in  the  bottom  layer  which  in  turn  depends  on  the  velocity  structure 
of  the  bottom.  To  obtain  a convenient  form  for  the  wave  function  ip,  the 
assumed  linear  velocity  profile  in  the  sediment  layer  is  replaced  by  a 
pseudollnear  one;  l.e.. 


o 


This  substitution  is  valid  since  most  of  the  energy  will  refract  within 
the  depths  where  the  pseudollnear  profile  is  a good  approximation  to  the 
linear  profile. 


. f 


The  use  of  the  pseudolinear  profile  means  that  the  wave  function 
in  the  sediment  layer  is  given  by  the  Airy  Function;  i.e., 

i|)(Z)  = Ai[C(Z)] 

In  this  case,  the  condition  specified  by  Eq.  IV. 20  implies 
<i»{0)  = Ai[iH0)]  = 0 


-(t)  ®)  ] = “i  • 


IV.  21 


where  is  a zero  of  the  Airy  F\mction. 

3.  Application  of  the  Method 

To  see  how  Hemna's  Method  is  applied,  consider  the  experimental 
transmission  loss  curve  pictured  in  Fig.  IV. U.  The  curve  was  produced 
from  data  taken  by  a bottomed  hydrophone  over  an  miconsolidated  sediment 
layer  about  1 km  thick. 

Fran  a knowledge  of  the  source  depth  and  the  water  profile, 
those  nulls  in  Fig.  IV. U iidiich  arise  because  of  Lloyd's  Mirror  Effects 
may  be  identified.  These  are  the  nulls  marked  A and  B in  Fig.  IV. U. 

The  niill  marked  E is  at  a range  where  the  direct  arrival  has  grazing 
incidence  and  hence  may  be  attributed  to  Eq.  IV. 19.  The  nulls  marked 
C €uad  D must  therefore  be  caused  by  refracted  arrivals.  Hanna  proposes 
a method  for  using  the  positions  of  these  nulls  to  evaluate  the  ratio 
C^/C  and  the  gradient  at  the  top  of  the  sediment  layer,  3. 
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FIGURE  1Z.4 

TRANSMISSION  LOSS  versus  RANGE 
(AFTER  HANNA) 
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From  the  range  coordinates  of  the  two  nulls  and  a knowledge 
of  the  water  profile,  the  grazing  angles  of  the  direct  arrivals  may  be 
determined.  With  these  grazing  angles  and  zeros  of  the  Airy  Function, 
one  has  two  equations  like  Eq.  IV. 21  which  may  be  solved  simultaneously 
for  C^/C  and  B.  When  Hanna  applied  this  technique  to  the  curve  of 
Fig.  IV. 4 he  obtained  the  following  results  which  were  in  good  agreement 
with  seismic  measurements , 


c /c  = 0.98 

o 


B = 1.5  sec' 


This  method  can  be  used  to  indicate  the  presence  of  velocity  gradients 
and  to  verify  or  confirm  results  obtained  by  other  methods.  In  situations 
in  which  propagation  loss  is  the  only  type  of  data,  this  method  might 
well  be  the  only  way  of  obtaining  anj'-  information  about  the  bottom. 
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APPENDIX  A 


DERIVATION  OF  EXPRESSIONS  FOR  INTERVAL  VELOCITY  AND  VERTICAL 
TRANSIT  TIME  IN  A LAYER  WITH  A CONSTANT  GRADIENT 


The  interval  or  layer  velocity  of  a layer  is  defined  to  be  the  time 

averaged  velocity  through  the  layer.  Although  a variety  of  different 

expressions  for  the  time  averaging  are  used,  the  definition  used  by 
13 

Clay  and  Rona  will  be  adopted  in  this  report.  Their  definition  is 
given  in  Eq.  A.l. 


The  quantity  AT  is  the  two-way  vertical  transit  time  in  the  layer,  and 
the  velocity  as  a function  two-way  time,  C(T),  is  referenced  to  the  top 
of  the  layer  in  question. 


Next,  consider  the  layer  to  have  a linear  sound  speed  variation,  i.e 


where  C^  is  the  sound  speed  at  the  top  of  the  layer,  g is  the  constant 
gradient,  and  Z is  the  vertical  coordinate  referred  to  the  top  of  the 
layer.  With  Eq.  A. 2,  the  vertical  two-way  time  to  the  depth  coordinate  Z 
may  be  calculated  as  follows. 


T = 2 


If  Eq.  A. 3 is  solved  for  Z,  one  obtains 


If  Eq.  A.^  is  substituted  into  Eq.  A. 2,  the  velocity  of  sound  as  a 
function  of  two-way  time  is  obtained,  i.e.. 


C(T)  - C e® 
o 


To  calculate  the  interval  velocity,  Eq.  A. 5 is  substituted  into  Eq.  A.l 
to  yield 

« 2 /-at 


-2-1 
AT  y 


'int  gAT 


Equations  A. 3 and  A. 6 are  used  extensively  in  the  text  of  this  repoirt. 


APPENDIX  B 


DERIVATION  OF  THE  RELATIONSHIP  BETWEEN  THE  RAY  PARAMETER 
AND  THE  DERIVATIVE  OF  THE  TRAVEL  TIME  CURVE 


A relationship  that  is  frequently  used  in  ray  theory  is  the 
following  one. 


The  quantity  p is  the  ray  parameter  of  a particular  path  defined  by 
Snell's  Law  to  be 


cos  9(Z) 
C(Z) 


The  derivative  in  Eq.  B.l  is  taken  along  a travel  time  curve  having 
fixed  source  and  receiver  depths.  It  must  be  remembered  that  this  is 
a derivative  along  a travel  time  curve,  not  along  a ray  path. 

Equation  B.l  may  be  obtained  by  considering  two  ray  paths  with 
infinitesimally  differing  ray  parameters.  This  config\iration  is  pictured 
in  Fig.  B.l.  From  the  geometry  of  this  configiiration,  in  the  limit  of 
infinitesimally  differing  paths,  one  may  write 

eo,  e(Z)  » 

dx 


COS  6(Z] 
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